Cosmic Microwave Background anisotropics with mixed isocurvature perturbations 
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Recently high quality data of the cosmic microwave background anisotropies have been published. 
In this work we study to which extent the cosmological parameters determined by using this data 
depend on assumptions about the initial conditions. We show that for generic initial conditions, not 
only the best fit values are very different but, and this is our main result, the allowed parameter 
range enlarges dramatically. 



PACS: 98.80-k, 98.80Hw, 98.80Cq 
I. INTRODUCTION 

The discovery of anisotropies in the cosmic microwave 
background (CMB) by the COBE satelhte in 1992 ^ has 
stimulated an enormous activity in this field which has 
culminated recently with the high precision data of the 
BOOMERanG [|, DASI || and MAXIMA-1 |] exper- 
iments. The CMB is developing into the most impor- 
tant observational tool to study the early Universe. So 
far, this data have however mainly been used to deter- 
mine cosmological parameters for a fixed model of initial 
fluctuations, namely scale invariant adiabatic perturba- 
tions |-|l|. 

In this work we investigate to which extent the deter- 
mination of cosmological parameters depends on initial 
conditions. Rather than exploring the most general pa- 
rameter space, we mainly want to show in a specific ex- 
ample how the allowed parameter range is enlarged when 
the requirement for purely adiabatic initial conditions is 
loosened. Therefore, in order to limit the computational 
effort we have chosen to vary some cosmological param- 
eters and keep the others fixed. We have set the to- 
tal density parameter f2tot = f^A + f^m = 1 and fixed 
r^A = 0.7 and fim = -I- ilb = 0.3. Here Ha denotes 
the density parameter due to a cosmological constant, 
f2A — A/S-Z/q, while f2c and flh are the density parame- 
ters of cold dark matter (CDM) and baryons respectively, 
and Hq = lOOh km s^^ Mpc~^ is the Hubble parameter 
today. For fixed Q\, Qc and spectral index ns = 1 we 
determine the parameters h and oj\, = Q^ih^ for generic 
initial conditions. 



II. INITIAL CONDITIONS 

A study considering the adiabatic together with just 
one isocurvature mode has been undertaken recently . 
To choose an even more generic set of initial conditions 
we follow the procedure outlined in Ref. ||l^ . The matter 
components present in the "standard" universe are CDM, 
baryons, massless neutrinos (v) and photons (7). Usually 
adiabatic initial conditions are assumed for the perturba- 
tions of these components. In Newtonian or flat-slicing 



gauge, this implies 

Ac = Ab = ^A^ = ^A^, and = Vh = = K, 

where A^ and Vi are the density fluctuation and the pecu- 
liar velocity potential of component i. However, this does 
not represent the most general set of possible initial con- 
ditions for cosmological perturbations. Apart from the 
adiabatic mode given above there are a baryon isocurva- 
ture mode (BI) a CDM isocurvature mode (CI) a neu- 
trino isocurvature density mode (NID) and a neutrino 
isocurvature velocity mode (NIV). The precise definition 
of these modes is given in Ref. |lj] and some observable 
consequences of deviations from a pure adiabatic model 
were already investigated in Ref. p^ . The most generic 
initial conditions are then given by a positive definite 
5x5 matrix M representing the amplitude of each of 
these modes, including all the possible cross-correlations 
(see Ref. ; a gauge-invariant formulation of the initial 
conditions is given in Ref. ||l^). We have noticed that 
implementing the initial conditions for all the modes in 
synchronous gauge was somewhat tricky, while this pro- 
cedure is very simple and numerically unproblcmatic in 
gauge-invariant formalism. 

For a fixed set of cosmological parameters we first 
compute the CMB anisotropy spectrum C]"* when only 
one of the elements of the correlation matrix is non-zero 
(Mij = 1, all other elements vanish) with a fixed spectral 
index ns = 1 for all modes. Next we set 



(1) 



2J = 1 



As already noticed in Ref. [|l7|], the BI and CI compo- 
nents of the correlation matrix are identical, up to an 
irrelevant multiplicative constant (which is due to the 
fact that fib 7^ ric). We have therefore restricted our 
analysis to the four modes AD, CI, NID, NIV without 
loss of generality. We then vary the correlation matrix 
M and the cosmological parameters o^b and h to search 
for the best fit to the data using a maximum likelihood 
method. 
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h = 0.65 




FIG. 1. CMB anisotropy spectrum for different values of 
the cosmological parameters LJh and h. We have shown the 
best-fit corresponding to a purely adiabatic case (dashed line) 
and allowing the most general set of initial conditions (solid 
line). The calibration and the beam size of the BOOMERariG 
data have been optimized to fit the mixed model (red/dark 
error bars) or the adiabatic model (green/light error bars). 
The parameter choice on top (tJb = 0.02, h — 0.65) can be 
fitted by both models while the values Uh = 0.042, h = 0.65, 
can only be fitted by a mixed model. 

Since the main point of this paper is to clarify the role 
of initial conditions and not to obtain the most realistic 
set of cosmological parameters, we restrict our analysis 



to the COBE § and BOOMERanG § data. For the 
BOOMERanG data we also take into account the calibra- 
tion and the beam size uncertainties |^ which we treat 
just like two additional (normally distributed) parame- 
ters of the problem. The best fits are computed using a 
downhill simplex method initiated after choosing a 
starting point randomly. Although this is certainly not 
the best method it proved to be sufficient for our pur- 
pose. In practice, each fit is the best fit among 15,000 
minimization runs with random initial conditions. The 
positive semi-definiteness of the correlation matrix M is 
insured by penalty functions on the sub-determinants of 
M (more details are given in jl^]). It turns out that the 
topology of the surface on our 14-dimensional param- 
eter space is quite complicated with many local minima 
and probably many degeneracies (see also the example 
discussed in |p^). 

In Fig. we show the best fit spectra for two different 
choices of the cosmological parameters Wb and h. Both of 
them are good fits if we allow for mixed initial conditions. 
On the plot we have also indicated the reduced x^. For 
a fixed choice of the parameters Wb, h the purely adia- 
batic model has only 3 parameters (the amplitude of the 
adiabatic mode, the BOOMERanG calibration and beam 
size). With 26 data points (7 from COBE and 19 from 
BOOMERanG) this leads to Fad = 26 - 3 = 23 degrees 
of freedom, while the mixed models have a symmetric 
4x4 matrix determining the initial amplitude, leading 
to a total of 12 parameters and hence only Fmix = 14 
degrees of freedom. If we also vary Wb and h, the number 
of degrees of freedom is lowered by 2. It is of course not 
surprising that for the parameters h = 0.65, ujh = 0.02, 
which are well fitted by the adiabatic model, the reduced 
of the adiabatic model is somewhat lower than the 
one of the mixed model. For the mixed model, the ab- 
solute is always lower than the one of the adiabatic 
model since the latter one is contained in the generic 
class of mixed models, but the reduced can of course 
be higher since it is divided by a smaller number of de- 
grees of freedom, Fmix < ^ad- This is actually the case 
for the parameters choice in the top panel of Fig. |^. 

To indicate what happens when models with mixed 
initial conditions are admitted we determine the likeli- 
hood functions of the cosmological parameters ujh and 
h by marginalizing over the initial conditions and the 
BOOMERanG calibration and beam size. The result 
is shown in Fig. || where the likelihood contours in the 
{ujb,h) plane for likelihoods of 50%, 68%, 95%, 99% are 
indicated for mixed models (top) and for purely adiabatic 
models (bottom). This figure is our main result. It is 
quite amazing to see to which extent the innermost good 
fit contour opens up once we allow for isocurvature com- 
ponents. Strangely, the only excluded region which re- 
mains is the upper left corner corresponding to the value 
of Wb inferred from big bang nucleosynthesis and the 
Hubble space telescope key project value for the Hubble 
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parameter (2^] of /i = 0.72 ± 0.08. 
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FIG. 2. The likelihood contours of 50%, 68%, 95%, 99% 
are indicated in the (ojb, h) plane for mixed models (top) and 
for purely adiabatic models (bottom) . The likelihoods are ob- 
tained by marginalization over the BOOMERanG calibration 
and beam size as well as over the initial conditions given by 
the matrix M for mixed models and by the amplitude of the 
adiabatic mode for adiabatic models. For the mixed model, 
the lowest corresponds to even higher values of ujh and h 
than those shown in the plot. 

There is absolutely no upper limit for Wb within the 
regime investigated here! Actually, this can be under- 
stood by noting that with adiabatic initial conditions, 
the strongest feature of a high baryon density universe is 
the asymmetry between even and odd acoustic peaks, but 
recalling that this asymmetry exists only in the matter 
dominated era and not in the radiation dominated era. 
With the parameters we have chosen here, the high value 
of 17 A implies a low matter content and hence, equiv- 
alence between radiation and matter happens not long 
before decoupling. For high values of h, the equivalence 
occurs earlier and the asymmetry of the peaks allows to 
put upper limits on the baryon density in purely adia- 
batic models. The second feature of high baryon density 



is to reduce the sound horizon and therefore to shift the 
acoustic peaks to the right. These two effects constrain 
the baryon density most strongly in the class of purely 
adiabatic models. In the case of mixed initial conditions, 
it is known that the presence of the isocurvature modes, 
especially the neutrino isocurvature modes, can not only 
shift the peaks significantly but also change their rel- 
ative heights (see lower panel of Fig. [l]). A high 
baryon density can therefore easily be accommodated in 
this framework. Note however that for high lu^ and low 
h, it is difficult to find a good fit, even for mixed models, 
because there is not enough power in the secondary peak 
region. The reason is that the early integrated Sachs- 
Wolfe effect boosts the first peak. 

We define the isocurvature content of a mixed model 
by a = (M22 + M33 + M44)/traceM, where Afn denotes 
the adiabatic mode amplitude. The isocurvature content 
in the good fit model shown in the top panel of Fig. |l| 
is only a — 0.12, while for the parameter choice in the 
bottom panel a — 0.69. Hence, if the cosmological pa- 
rameters are close to those chosen in the top panel, we 
can conclude that the cosmic perturbations are predom- 
inantly adiabatic. In Fig. ^ we show the isocurvature 
content a of the best fit model obtained by minimizing 
by variation of the initial conditions at fixed values 
of the cosmological parameters. Clearly, the further we 
move away from the parameter region well fitted by the 
purely adiabatic model, the higher becomes the isocur- 
vature contribution needed to fit the data. 
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FIG. 3. The isocurvature content a of the best fit mixed 
model at fixed parameter values {cuhjh) is indicated. The 
contours a = 0.2 to 0.9 in steps of 0.1 are shown. 

The main non-adiabatic component of our best fits is 
the NID mode. This was to be expected as this mode and 
its correlator with the adiabatic mode can shift the peak 
positions and can substantially add or subtract from the 
second peak. It is known that for interacting species the 
non adiabatic part of the perturbations tends to decay 
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with time. Therefore the generation of this NID com- 
ponent can only occur after neutrino decouphng, that 
is at T < 1 MeV. Whether or not such a phenomenon 
can occur at low energy is an open question. However, a 
neutrino isocurvature perturbation can also be due to a 
fourth species of sterile neutrinos which may have decou- 
pled very early in the history of the Universe. The same 
remark applies of course also to the CDM isocurvature 
mode. Note that the energy density of this fourth neu- 
trino type cannot be very high in order not to contradict 
the light element abundances, but there is nothing which 
prevents (at least in principle) the presence of large per- 
turbations in this fluid. 

III. CONCLUSION 



dipole Vy rather than the photon density perturbation 
A^, these two quantities depending in a different way on 
the initial conditions. In the same vein, using the normal- 
ization of the matter power spectrum (provided it can be 
measured accurately...) also helps to break some of the 
degeneracies induced by the isocurvature modes. 
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We have shown that allowing for isocurvature contri- 
butions one may very well fit present CMB data with 
a set of cosmological parameters which differs consider- 
ably from the one preferred by adiabatic perturbations. 
More important, allowing for generic initial conditions, 
the ranges of cosmological parameters which can fit the 
CMB anisotropy data widen up to an extent to become 
nearly meaningless. 

On the other hand, assuming measurements of cosmo- 
logical parameters from other methods like direct mea- 
surements of the Hubble parameter which yield h ~ 0.65 
and big bang nucleosynthesis which implies Wb ^ 0.02, 
we can use the CMB to limit the isocurvature contribu- 
tion in the initial conditions (or other unconventional fea- 
tures) and thereby learn something about the very early 
universe, the inflationary phase which probably has gen- 
erated these initial conditions. 

We believe that it is already interesting to note that for 
cosmological parameters in the range preferred by other, 
CMB independent, measurements (JIa ~ 0.7, fim ^ 0.3, 
h 0.65, ujh ~ 0.02) the isocurvature contribution in the 
initial conditions has to be relatively modest. 

Finally, our work shows the danger of calling parameter 
estimation by CMB anisotropy experiments a "parameter 
measurement" since the results depend so sensitively on 
the underlying model assumptions. We rather consider 
CMB anisotropics as an excellent tool to test model as- 
sumptions or consistency. In the light of these findings, 
the importance of non-CMB measurements of cosmolog- 
ical parameters can clearly not be overstated. In short, 
CMB seems to be a better tool to investigate the pri- 
mordial parameters (i.e., the initial conditions) rather 
than the cosmological parameters (fi's, h etc). Hope- 
fully, measuring the CMB polarization will represent an 
additional non-trivial mean to remove the degeneracy 
present between cosmological parameters and initial con- 
ditions. This has been studied in detail in Ref. ijl^. The 
main reason for this is that polarization, which is induced 
only during decoupling, is mostly sensitive to the photon 
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